ENDPOINT ESTIMATES AND GLOBAL EXISTENCE FOR THE 
NONLINEAR DIRAC EQUATION WITH POTENTIAL 
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Abstract. We prove endpoint estimates with angular regularity for the wave 
and Dirac equations perturbed with a small potential. The estimates are 
applied to prove global existence for the cubic Dirac equation perturbed with 
a small potential, for small initial data with additional angular regularity. 
This implies in particular global existence in the critical energy space for 
small radial data. 



1. Introduction 

The main topic of this paper is the cubic massless Dirac equation on 
turbed with a potential 
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iut — + Vix)u + Psiu, u), m(0, x) = f{x), 
p3 , 



(1.1) 
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and P3{u,u) is any homogeneous, C -valued cubic 



polynomial. We recall that the Dirac operator V is defined as 
where 9, are the partial derivatives on Ri^ and aj are the Dirac matrices 
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(1.2) 



The anticommutation relations 

aiau + auai = 2Skil4 

imply that = — /4A is a diagonal operator, showing the intimate connection of 
the massless Dirac system with the wave equation. 
The unperturbed nonlinear Dirac equation 



iut ~ Vu + F{u), u(0, x) — f{x) 



(1.3) 



is important in relativistic quantum mechanics, and was studied in a number of 
works (see e.g. [20], [7], [19], [18], [8], [17], [IG] and for the more general Dirac- 
Klein-Gordon system see [G], [5]). In particular, it is well known that the cubic 
nonlinearity is critical for solvability in the energy space H^] global existence in 
is still an open problem even for small initial data, while the case of subcritical 
spaces i?", s > 1 was settled in the positive in [8], [17]. 

Criticality is better appreciated in terms of Strichartz estimates, which are the 
main tool in the study of nonlinear dispersive equations. The identity 

e^*^/ = cos{AD\)f + ^^^^^P/, \D\ = (-A)i/2 



Date: January 15, 2013. 



1 



NONLINEAR DIRAC EQUATION 2 

shows that the estimates for the Dirac flow e''^ are immediate consequences of the 
corresponding estimates for the wave flow e'*l^l, with the same indices, restricted 
to the special case of dimension n = 3. For the wave equation on R^"*"", n > 3, 
Strichartz estimates can be combined with Sobolev embedding and take the general 
form 

|||i5|-+^-fe^tPI/||^,^.< 11/11^, (1.4) 

for all p, r such that 

1 1 2 

pe[2,oo], 0<-< 



r 2 {n — l)p 
We are using here the mixed time-space L^L^ norms defined by 

\\u{t,x)\\LPLi = ||||u|L«|Lp • 

Notice that the limiting case r = oo 

ie^'I'^l/llL^L- <III^|'^/I|l^ (1.5) 



is always excluded and is indeed false for general data. See [If] and [f3] for the 
general Strichartz estimates; concerning the limiting case r = oo, see [f4], [9]. The 
corresponding estimates for the Dirac equation are given in [3]. 
In particular, the endpoint estimate 

l|e"^/||L-^L~<|||i?|/IU^ (1.6) 

fails. To see the connection with the critical equation (1.3), we rewrite it as a fixed 
point problem for the map 

v^^{v) = e''^f + t f e'(*-*')^P3(w(i'))di'. 
Jo 



If (1.6) were true one could write 







+00 



and in conjuction with the conservation of energy this would imply 

\\^{y)\\LrHi + ||$(w)||l?l- < WfWm + \\v\\L^m\\v\\l2L'>^. 

In other words, a contraction argument in the norm || • 11^2^00 + || • ||^oo^i would be 
enough to prove global existence of small solutions to (1.3). 

It was already noted in [14] that (1.5) is true for radial data when n — 3. 
This remark is not of immediate application for the Dirac equation, since solutions 
corresponding to radial data need not be radial (due to the fact that the operator 
V does not commute with rotations of M'^). Nevertheless, for radial and even 
more general data, in [16] global existence was achieved via finer estimates, which 
separate radial from angular regularity. We introduce the natural notations 

and 

ll/lli-L^ =SUp||/(r • )||Li.(s„-i). 

Then the following estimate is proved in [IG]: 

^ = 3, ||e"l^l/|U2,^i£ < v^.|llDl/lU2, Vp<(X3. (1.7) 
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This gives a bound for the standard L^L°° norm via Sobolev embedding on the 
unit sphere 



Jt\D\ 



/II 



< l|ALe**l^l/|| 



< 



l^|AL/ll 



L2, 



2 

P> - 
e 



(1.8) 



where the angular derivative operator A* is defined in terms of the Laplace-Behrami 
operator on §"~^ as 

Af, = (l-As„-i)^/2_ 

Using (1.8) one can prove global existence for (1.3) provided the norm || |Z)|A5/||l2 
of the data is small enough for some s > 0. In particular, this includes all radial 
data with a small norm. 

Our main goal here is to extend this group of results to the equation (1.1) 
perturbed with a small potential V{x). We consider first the linear equation 

iut =Vu + V{x)u + F{t,x). (1.9) 

The perturbative term Vu can not be handled using the inhomogeneous version of 
(1.7) because of the loss of derivatives. Instead, we prove new mixed Strichartz- 
smoothing estimates (Theorem 2.3) 
t 

e^(*-^)l^li^(s,a;)ds <\\(x)i+\D\ 



n>3 
where 



LiLl 



for n — 3, 
for rt > 4, 







n 

a = 1 . 

2 



(1.10) 



(1.11) 



Remark 1.1. As a byproduct of our proof, we obtain the following endpoint esti- 
mates for the wave flow with gain of angular regularity (Theorem 2.1): 

n>3, ||e'*l^l/L?L^LS < IIA^/II^:^ (1.12) 

where a is as in (1.11). Although this was not the main purpose of the paper, it 
is interesting to compare (1.12) with known results. In dimension n = 3, estimate 
(1.12) is just a special case of Theorem l.l-III in [IG] where (1.12) is proved with 
cr = — |; it is not known if this value is sharp, however in the same paper it is 
proved that the estimate is false for cr < — |. On the other hand, to our knowledge, 
estimate (1-12) for n > 4 and (1-10) for n > 3 are new. The literature on these 
kind of estimates is extensive and we refer to [10], [12] and the references therein 
for further information. 

Combining (1.10) with the techniques of [4] we obtain the following endpoint 
result for a 3D linear wave equation with singular potential. Analogous estimates 
can be proved for higher dimensions; here we chose to focus on the 3D case since 
the assumptions on V take a particular simple form: 

Theorem 1.1. Let n = 3 and consider the Cauchy problem for the wave equation 
utt - Au + V{x)u = F, u{0,x) = f{x), ut{0,x)^g{x) 

under the assumptions: 

(i) V{x) is real valued and the positive and negative parts V± satisfy 

C S 



V- < 



\x\2-' + \x\^ \x\^-' 

for some S, e sufficiently small and some C > 0; 
(ii) —A + V is selfadjoint; 



(1.13) 
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(iii) is not a resonance for —A + y_ (in the following sense: if f is such that 

(-A + V-)f = and (x)-^ f e L^, then f = 0). 
Then the solution u(t, x) satisfies the endpoint Strichartz estimate 

Ml-^l^li < II/IIhi + llffllL^ + \\{x)i+F\\L2L.. (1.14) 

The next step is to prove suitable smoothing estimates for the Dirac equation 
with potential 

iut ^'Du + V{x)u + F{t, x) 

(see Proposition 4.2 and Corollary 4.3). Then by a perturbative argument we obtain 
the following endpoint estimates for the linear flows: 

Theorem 1.2. Assume that the hermitian matrix V{x) satisfies, for S sufficiently 
small, C arbitrary and a > 1, with Va-{x) ~ 1 log + (x)^^"' , 

\V{x)\ < \VVix)\ < (1.15) 

Then the perturbed Dirac flow satisfies the endpoint Strichartz estimate 

\\e''^''^''^f\\L^LrLl<\\f\\m- (1.16) 
If the potential satisfies the stronger assumptions: for some s > I, 

\\KV{\x\ ■ )|U.(s.) < \\K^V{\x\ ■ )|U.(s.) < (1.17) 

then we have the endpoint estimate with angular regularity 

l|A:e^*(^+^)/llL?L;^Li < WKfWm (1.18) 
and the energy estimate with angular regularity 

\\Ke'''-'^^''^f\\L-m<\\Kf\\m (1.19) 

We can finally apply Theorem 1.2 to the nonlinear equation (1.1) and we obtain: 

Theorem 1.3. Consider the perturbed Dirac system (1.1), where the 4x4 matrix 
valued potential V{x) is hermitian and satisfies assumptions (1.17). Let Psiu^u) 
be a -valued homogeneous cubic polynomial. Then for any s > 1 there exists Cq 
such that for all initial data satisfying 

\\Kf\\m<eo (1.20) 
the Cauchy problem (.5.1) admits a unique global solution u G CH^ H L^L°° with 

In particular, problem (1.1) has a global unique solution for all radial data with 
sufficiently small norm. 

Remark 1.2. It is clear that our methods can also be applied to nonlinear wave 
equations perturbed with potentials, and allow to prove global well posedness for 
some types of critical nonlinearities. This problem will be the object of a further 
note. 

Remark 1.3. We did not strive for the sharpest condition on the potential V , which 
can be improved at the price of additional technicalities which we prefer to skip 
here. Moreover, the result can be extended to more general cubic nonlinearities 
\F{u)\^\u\\ 

Notice also that we need an angular regularity s > 1 on the data, higher than 
the s > assumed in the result of [Ki]. It is possible to relax our assumptions to 
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s > 0; the only additional tool we would need to prove is a Moser-type product 
estimate 

\\K{nv)\\Li < WuU^WKvUi + ||A>|U. s > 

and an analogous one for Af^\D\{uv). This would require a fair amount of calculus 
on the sphere and here we preferred to use the conceptually much simpler 
algebra property of H^{S"~^) for s > ^^-j-'-- 

On the other hand, the extension of our results to the massive case 

iut = Vu + V{x)u + mf3u + F{u), m ^ 

requires a different approach and will be the object of further work. 

Acknowledgments. We are indebted with Ilia Krasikov and Jim Wright for invalu- 
able conversations which helped shape up the proofs in Section 2. 

2. ENDPOINT ESTIMATES FOR THE FREE FLOWS 

To fix our notations, we recall some basic facts on spherical harmonics (see [22]) 
on M", n > 2. For A: > 0, we denote by Hk the space of harmonic polynomials 
homogeneous of degree k, restricted to the unit sphere S"~^. The dimension of Hk 
for fc > 2 is 

while do = 1 and di — n. T-Lk is called the space of spherical harmonics of degree 
k, and we denote by F^', I < I < dk an orthonormal basis. Since 



k=0 

every function f{x) ~ f{ruj), r ~ \x\, can be expanded as 



/w = EE/feW^feH (2.1) 

k=0 1=1 

and we have 

fc>0 

l<i<d;, 

where we use the notation = L^(S^~^). More generally, if Ag is the Laplace- 
Beltrami operator on S'^~^ and 

A„ = (l-As)i/2, 

we have the equivalence 

k>0 
l<l<df. 

As a consequence we have the equivalence 

l<l<df^ 

In a similar way 

l|V/||i.(R„) = (-A/,/)i. 

- E (lk"^5^/^(Olli?(o,oo)+fc'lk"^'/^(-)llL.(o.oo)) (2-3) 
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where we used the foUowing representation of the action of A 



More generaUy we have for integer m 
-A(l~As)'"/(^) = ^(l+fc(fc+n-2))"%' 

which imphes 

l|VA-/||i.(j,„, = (-A(l-A5r/,/U. 



Y: (fc)^'" {h^drfim.^o.oo, + fc^!|r^/^(r)|li.(o,oo)) (2.4) 



k>0 
KKdf, 



and by interpolation and duahty we see that (2.4) holds for all m e 
We shall estimate the solution using the following norm: 

II/IIl-lJ = sup ||/(rw)||i2 (s„-i). 

r>0 

Theorem 2.1. For all n> A the following estimate holds: 



711 



while for n = 3 we have 



/II 



< 



(2.5) 



(2.6) 



Remark 2.1. In dimension n = 3 the previous result is a special case of the stronger 
estimate proved in [16]: 



||e"l^l/|| 



< 



(2.7) 



Notice that it is not known if estimate (2.7) is sharp. For higher dimension, estimate 
(2.5) seems to be new; it is reasonable to guess that this result is not sharp and 
might be improved at least to 



/II 



LjLS^L 



2 < \\K 



711^^1 



e > 0. 



(2.8) 



Proof. It is well known that the Hk spaces are invariant for the Fourier transform 
J^, and more precisely 



j-(c(r)y,'(c.)) iO^giimk 

where g is given by the Hankel transform 



g{r) = (27r)^i ''r ' 



c{p)J^^r^{rp)p'idp. 



(2.9) 



(2.10) 



Here Ji, is the Bessel function of order which we shall represent using the Lommcl 
integral form 



My) 



7r5r(^+ 1/2) 



(2.11) 



Now, given a function f{x), we denote by / its inverse Fourier transform and with 



fi^{r) the coefficients of the expansion in spherical harmonics of /: 

oo dk 

/ = EE/fcW^'^'H- 

fe=0 1=1 



(2.12) 
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Recalling (2.9) we obtain the representation 



/(x)=5](27r)tr^-|x|i-tr,'(^^)^°°/^(p)J,+ ^(Np)ptdp (2.13) 



which implies 

e'*l^l/ = ^(2^)t 



jtp fi 



fUp)Jk+^i\x\p)p^dp. (2.14) 



Consider now Lommel's formula (2.11) for Ji,; since e'^^ = (iy) '^d^{e'^^y), after k 
integration by parts we obtain 



J, 



{y)^cuy^-' j ((1 - AY^"^) 



dX 



with 



Thus we can write 

/•OO 



Cfc 



7r5r(S^ + fc)' 



(2.15) 



(2.16) 



e'"'fi{p)Jk+^^{\Ap)p^dp 







c, / ^an(i-A^) 



l+(p)/^(p)p""'e^''^*+^l-l)dp 



(2.17) 



where l+(p) is the characteristic function of (0,+oo); regarding the inner integral 
as a Fourier transform we arrive at 



Ck dl ((1 - gi{t + X\x\)d\ 



(2.18) 



where 

gi{p) = l+{p)fiip)p^^-' 
In conclusion, we have the following representation 

e^t\D\j^Y^^2^yi^-kyi(^_^^^^ |'^9^-((l-A2)'^+^)?i(t + AN)dA (2.19) 

where the constants Cfc are given by (2.16) and gl. by (2.18). Notice that similar 
representations play a fundamental role also in [10], [12] In particular this gives for 
the norm of the solution at i, |a::| fixed the formula 



iie^*i^i/(ix|.)iii. c^Ei^^-i' y /'((i-^')'"*-"^ 

We now need the following estimate: 
Lemma 2.2. Let Qk{x) be the function 

((l-a;2)'=+^) 



gi{t + X\x\)dX 



(2.20) 



Qkix) = 

Then we have on x ^ 



\Qkix)\ < (fc)i-* z/n>4, \Qk{x)\ < 1 z/n = 3. 



(2.21) 
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Proof. We recall that the Jacobi polynomials are defined by 

pi"'^(-) = - + -r'^ [(1 - -r^'ii + . (2.22) 

We shall use some standard properties of these polynomial which can be found in 

k 



[1]. The function Qk can be expressed in terms of P\°''°'\x) with a — (n — 3)/2 as 



fc!(l-a;2)'^ 



1-3 n-3\ 



Thus in order to estimate Qk we need a bound for the function 

Ux) = {l-x'rPi^^^\x), a=^. (2.24) 

The following approach was suggested by Ilia Krasikov, see [15]. Consider the 
second orderd differential equation 

r{x)+p{x)r{x)+q{x)f{x)^0 

on the interval (—1, 1), and define the Sonine function as 

s(f,x)=f{xr + tif)l 

under the assumption q > 0. It is easy to check that function S satisfies the relation 

The function Ta{x) defined in (2.24) satisfies the differential equation 
, 2(2a-l) , (fc + l)(2a + fc) , , 

1 — a;^ 1 — x'^ 

so that the associated Sonine function 

Sa{x) = Tl + (fc^i)(2a + fc)^'^' 

satisfies 

From this identity it is clear that Sa has a maximum at a; = provided a > 1/2 
i.e. 71 > 4. In this case we have 



Sa{x) < SaiO) = TaiOf + , , = Pi^'''\0)' 



(fc + l)(2a + /c) ^ ^ (fc + l)(2a + fc) 

Now we recall that, for even k > 2, 

P^'^)(0) = r ^ {-l)''k-^ P^")'(0) = 

^ ^ (_2)fcr(| + i)r(| + a + i) ^ ' ^ ' 

where we used the Stirling asymptotics 

fc! ~ fc'=^i/2g-fc^ r(fc + a + 1) ~ fc'^+^'-i/^e-'^. 
In a similar way, for odd fc, 

pia.a) _ Q pia.aY.Q. _ T (fc + a + 1) (-1)^-1 

(U) - u, (U) - (_2),_ip (I + 1) r (I + - n - ^ 

Thus for all values of fc > 1 we have 

1 



|T,(x)| < v/5:(^< — 
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and by (2.23) we conclude that, for fc > 1 and |a:;| < 1, 

\Qk{x)\<k^-'i 

which is precisely (2.21) for n > 4. 

In the remaining case n — S we have a — and the best we can do is to use the 
sharp inequality IP^^'""*! < 1 to obtain 



iQ.(-)i=|iipri<i- 



□ 



Using the Lemma, we can continue estimate (2.20) as follows 



l|e"i^i/(N-)llii<E-' {l\ 



giit + X\x\)\dX 



where 

ujk^l ifn = 3, tjfe = (fc)i~* ifn>4. (2.26) 

Since 



\t{t + X)\dX<M{gi)it) 



l'jg^^{t + X\x\)\dX^^J'^ 

where M(g) is the centered maximal function, we obtain 

l|e^*l^l/(N.)||i. <E-^^(5i)(ir- 
Now we can take the sup in which gives 

||e^*l^l/|li.^. <E..^M(5l)(^)^ 
and integrating in time, by the boundedness of the maximal funcion, we obtain 

It is immediate to check that the last sum is equivalent to 

where a — 1 —n/2 for n > A, which proves (2.5), and a — for n = 3, which proves 
(2.6). □ 

Although the method of proof of Theorem 2.1 is probably not sharp for the 
homogeneous operator, it has the advantage that it can be adapted to handle 
also the nonhomogeneous term and gives the following mixed Strichartz-smoothing 
estimate: 

Theorem 2.3. For any n > 3, the following estimate holds: 



3'(*-^)l^lF(s,x)ds 



< ||(x)H|i?|^ASFL.i. (2.27) 



where 

Ti 

(7 = 1 — — if n > A, a ~ if n = 3. 

Proof. As in the proof of the previous theorem, we expand F in spherical harmonics 
and we obtain the representation 

/ e^(*--)l^lF(s,x)ds = 
Jq 

= J2i2n)U-'^c,Yl (^^) J'^ d', ((1 - A^)'=+^) Gi{s,t- s + X\x\)dX 

(2.28) 
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with the constants Ck as in (2.16), where the functions G^, are defined as follows: 
denoting by Fl{t,r) the coefficients of the expansion into spherical harmonics of 
the inverse Fourier transform F = J-^^{F) 



and by Gj, the functions 



the r) are the Fourier transforms of G^. in the second variable: 

/'+00 



Thus applying Lemma 2.2 we obtain 



fI J-\x\ Jo 



(2.29) 



t-s + X)\ (2.30) 



where Wk is the same as in (2.26). We estimate the integral in s as follows 



\G[\ds< 



+00 



(A + t - s)5+(A + t - s)-^-\G[{s, \ + t-s)\ds 



< 



{X + t- sy+\G[{s, X + t- s)\'ds = Ql.(A + t), 



where 



\Gi{s,t,-sr{fi~sy+d.s 



Thus we see that 

r\x\ rt 



- r dX r ds\GUs,t-.s + X)\<^ f^\i{X + t)dX<M{Qi){t) 

^\ J~\x\ Jo Fl J-\x\ 



Coming back to (2.30) we obtain 



5'(*-'*'l^lF(s,a;)ds 



and taking first the norm, then the sup in \x\, then the Lj norm, by the 
boundedness of the maxiaml function we have 

2 



e'(*-'^)l^li^(s,a;)ds 



< 



J2^i\\Qim 



The definition of Qj. implies 



\Qiit)\^dt = J J \Giis,^,-s)\■'{^l~sy+dsd^,^\\Giit,r){r)^- 
and hence 



2 



2*(*-^)l^lF(s,:E)di 



< 



LiL^L^, 
I 



E 



^fcl|Gl(t,r)(r)H||2 



(2.31) 



Recalling the definition (2.29) of Gj,, we see that to obtain (2.27) it is sufficient to 
prove the following general inequality for s = 1/2+ and arbitrary a: 

E(fc>'^ ||(y>^-^A^. {i+{x)y'-'fUm% ^ Wi^rKlDl'^fWhm- (2.32) 

Here as usual fl denotee the coefficients in the expansion in spherical harmonics of 
the inverse Fourier transform / = T^^ j . 
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First of all, since T ^d-D] ^ /) = |^| 2 j ^ we see that it is enough to prove, for 
< s < 1 and arbitrary cr, the slightly simpler 



< 



LI 



(^)'*^w/llL2(Ri')- 



(2.33) 



The inequality will follow by interpolation between the cases s = and s = 1; 
indeed, we can regard it as the statement that the operator T defined as 



T 



: / ^ [{yYFx^y (l+(A)A^.gl.) }^ ^ , g = f 



which associates to the function / the sequence of coefficients in the expansion of 
J""i(A-'"/), multiplied by X'-"'~^^/^l+, transformed again and multiplied by {y}^, 
is bounded between the weighted spaces 



When s = we have by Plancherel's Theorem and by (2.2) 

2 



Y,{kr^\\T,^y (l+(A)A"^/^(A)) 



J2{kr^ A^/^(A) 



\Kf\\ 



ii(0,A) 

Since A;^ commutes with the Fourier transform, indeed 

again by Plancherel we obtain (2.33) for s = 0. 

To handle the case s = 1 we consider the quantity 



yTx^y[l+{\)\-^Jl{\) 



LI 



< 



Lj(0,oo) 



1+(A)A^/^(A) 



L5(0,oo) 



< 



A^/^(A) 



Multiplying by (fc)^"^, summing over l,k and recalling (2.4), we obtain 



Y^ikf^ yT,^y(l+{X)X'^fiiX)) ^^<\\VAZf\\h 



ii(0,oo) 



where the last term can not estimated by (2.4) because of the factor fc^ which 
vanishes when fc = 0. However we have 



/o°(A) 



f{XLu)dX 



1^1=1 



AZf{XLo)dX 



which implies, using Hardy's inequlity 



A"^/o°(A) 
Thus we have proved 



< 



L|(0,oo) 



Kf 



1^1 



< l|VA^/|U.. 



L2(R") 



J^ikf" y^x^y{l+iX)X^flw) ^,<\\yKf\\h 



\x\Kf\\ 



L2 



again by the commutation of A;^ with the Fourier transform. This gives (2.33) for 
s — 1 and concludes the proof of the Theorem. □ 
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Remark 2.2. Since the operator A^^ commutes with \D\, estimates (2.5), (2.6) and 
(2.27) obviously generahze to the fohowing; for any real s > 0, 



and 



where 



F(s, x)ds 



<\\(x)-- + \D\ — K+''F\ 



(2.34) 



(2.35) 



(7=1 



if n > 4, 



(7 = if n = 3. 



From the previous estimate for the free wave equation it is not difficult to obtain 
analogous endpoint Strichartz and Strichartz-smoothing estimates for the 3D Dirac 
system: 

Corollary 2.4. Let n = 3. Then the flow e**^ satisfies, for all s >{), the estimates 



and 



Af 



s^(*-*')^F(i',x)di' 



< 



Proof. If u solves the problem 

iut+Vu = Q, it(0) = /(x), 
by applying the operator {idt — 2?) , we see that u solves also 

□u = 0, u(0) = /(x), ut{Q) = iVf. 
This gives the representation 



itT) 



f = cosit\D\)f + i 



sin(t|D|) 
\D\ 



Vf. 



(2.36) 
(2.37) 

(2.38) 
(2.39) 

(2.40) 



Moreover, we recall that the Riesz operators are bounded on weighted 

spaces with weight (x)" for a < n/2. Thus in the case s = estimates (2.36), (2.37) 
are immediate consequences of the corresponding estimates for the wave equation 
proved above. 

In order to complete the proof in the case s > 0, we need analyze the structure 
of the Dirac operator V in greater detail. Following [23], we know that the space 
L^(M'^)^ is isomorphic to an orthogonal direct sum 



L2(]R3)4 



i'(O,+oo;dr)0iJ,„^,fe^.. 



1 = 1 s. 1 

J 2 ' 2 



±(3 + 1/2) 

Each space Hmj.kj has dimension two and is generated by the orthonormal basis 
{^m fc ' fe }' which can be explicitly written in terms of spherical harmonics: 
when fcj = J + 1/2 we have 



V2j 
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/27 



V 



while when kj = — + 1/2) we have 



m j , kj 



.,-1/2 



V 








mj+1/2 



The isomorphism is expressed by the exphcit expansion 



/•CSO 





]dr. 



(2.41) 
(2.42) 
(2.43) 

Each i^(0, +oo; dr)®Hm.^kj is an eigenspace of the Dirac operator I? = ^ ajSj 
and the action of T) can be written, in terms of the expansion (2.41), as 



with 



Notice also that 



2 I ^nij 



\, dr "^^'''^ r ^'"j'^^j / r \dr " 



From decomposition (2.41) it is clear that the operator A^, which acts on spherical 
harmonics as 

azy;^ ^ {1 + £{£ + ■ y;^ , (2.44) 

does not commute with V. Indeed, each space Hmj,kj involves two spherical har- 
monics Y^ with two values of £ which differ by 1, and T) swaps them. However, 
the modified operator defined by 

obviously commutes with P, thus estimates (2.36), (2.37) are trivially true if we 
replace A with A. It remains to show that we obtain equivalent norms. The 
equivalence 

WJWli ^ WJWli 

follows directly from (2.44), (2.4.5) and (2.43). Moreover, A and A commute with 
A, hence with and this implies 

|||i^|A^/|U.^|||i?|A^/|U. 

or, equivalently, 

i!a:,/ii«.^iia:,/ii^,. 

This is sufficient to prove (2.36). Since A and A also commute with radial weights 
we have 



II(^)'+I^|a:/|| 



L2 



L2 



which gives (2.37). 



□ 
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3. The wave equations with potential 

Our next goal is to extend the results of previous section to the case of per- 
turbed flows. This will be obtainfed by a perturbative argument, relying on the 
smoothing estimates of [4] and the mixed Strichartz-smoothing estimates of the 
previous section. In [4] smoothing estimates were proved for several classes of dis- 
persive equations perturbed with electromagnetic potentials (while the ID case was 
analyzed in [2]). For the wave equation in dimension n > 3 the estimates are the 
following: 

Proposition 3.1. Let n > 3. Assume the operator 

-A + W{x, £>) = -A + a{x) ■ \7 + bi{x) + 62(2;) 
is selfadjoint and its coefficients satisfy 

- Ni- + |.p|iogN|<^ 

for some S,e > sufficiently small and some a > 1/2, C > 0. Moreover assume 
that is not a resonance for — A-I-&2- Then the following smoothing estimate holds: 

+ |x|)-ie**^^^/|U.i. < ll/IU.. (3.3) 

The assumption that is not a resonance for — A -I- 62(2^) here means: if (—A + 
62)/ = and {x)'^f e L2 then / = 0. 

Combining Proposition 3.1 with (2.27) we obtain the following Strichartz end- 
point estimate for the 3D wave equation perturbed with an electric potential: 

Theorem 3.2. Let n — 3 and consider the Cauchy problem for the wave equation 

utt - Au + V{x)u = F, w(0,x) = /(x), UtiO,x)=g{x) 

under the assumptions: 

(i) V{x) is real valued and the positive and negative parts V± satisfy 

V+ < , V- < (3.4) 

" \x\i-' + \x\^ ~ |a;|^-^ + |a;|2 

for some S, e sufficiently small and some C > 0; 

(ii) —/S. + V is selfadjoint; 

(iii) is not a resonance for — A + 

Then the solution u(t, x) satisfies the endpoint Strichartz estimate 

M^L^Li < m + llfflU^ + 11(2;)^+^^.^.. (3.5) 
Proof. We represent u{t,x) in the form 

u{t,x) =L + IL -LLI 



where 



and 



/ = cos(t|D|)/ + ^^5, 

LL= [ \D\-^ sm{{t - s)\D\)F ds, 
Jo 

111= I \D\^^siTi{{t- s)\D\)Vuds. 
Jo 
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We can use (2.5) to estimate / and (2.27) to estimate // in the norm L^L^L^ 
directly. On the other hand, applying (2.27) to /// we get 

WlIIh-L^Ll < \\{x)-^ + Vuh2L-2 < \\{x)-^+T,V\\L^\\T-'uh2L2. 

By assumption {x)^^t^V is bounded on M", moreover we are allowed to use (3.3) 
since V satisfies the assumptions of Proposition 3.1. Notice that (3.3) implies 

iirr'"iu^L^ <ii/iiL^ + iii^r\9iiL^ 

and in conclusion we have proved 

Wlllh^^LrLi < \\f\\L-^ + \\\D\-'g\\L'^ 

which completes the proof of (3.5). □ 

Analogous estimates can be proved for the Klein-Gordon equation, or in higher 
dimension n > 3, for first order perturbations, and for angular derivatives of the 
solutions. We omit the details since we prefer to focus on the Dirac equation here. 

4. The Dirac equation with potential 

We consider now the perturbed Dirac operator T) + V{x) where V{x) is a small 
4x4 hermitian matrix valued potential. We prove here more general versions of 
the estimates given in [4] , [3] in order to include angular regularity. We begin with 
the free Dirac equation: 

Proposition 4.1. The free Dirac flow satisfies, for all a > 1 and s > 0, the 

smoothing estimates (with Wa{x) — \x\{l + | log jxjD'^J 

Ww-'/'Ke^'-^fh.^. < WKfW^. (4.1) 

and 

-1/2^3 f\Kt-t')Vp^^.^^^> <\\n^l/^At,F\\^.^. (4.2) 

Proof. When s = 0, both estimates follow from the resolvent estimate 

||w-i/2i?p(z)/|U.(R3) < CWw'J^fh^iM^), z ^ M, 

with a constant uniform in z, proved in [3] using a standard application of Kato's 
theory (see also [4]). The case s > is proved exactly as in Corollary 2.4, first 
by replacing A^^ with A^ which commutes with the flow, and then by using the 
equivalence of norms. □ 

We consider now the case of a perturbed Dirac system 

iut — Vu + Vu 

where V{x) is a 4 x 4 matrix potential. If V is hermitian and its weak norm 
is small enough, the operator V + V is selfadjoint as proved in [3]. In all of the 
following results the assumptions on the potential are somewhat stronger than this, 
so in all cases the unitary flow e***^^+^^ will be well defined and continuous on 
L'^{R^y by spectral theory. 



Proposition 4.2. Let V{x) be a hermitian 4x4 matrix on such that 

\V{x)\<^^, w.{r)=r-{l + \\ogr\r (4.3) 

Wa{\x\) 

for some S > sufficiently small and some ct > 1. Then the perturbed Dirac flow 
^itiv+v) gfj^iigjigg ^/jg smoothing estimates 
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-1/2 j e^{t-t'){■D+v)p^^^r■^^^, 

" some s > 1 the 



<l 

If in addition V satisfies for some s > 1 the condition 

S 



LfLl 



Wa{r) ' 

then we have, for all < s <2, the estimates with angular regularity 

-1/2^. f\^(t-t'Kv+v)F^tr^dt' < lU^'KFW,.,.. 
Jo 

Proof. If u solves 
we can write 



(4.5) 

(4.6) 

(4.7) 
(4.8) 



iut ^Vu + Vu + F, u{0) = / 



Jtv _ 



f + i e'^'-'">'^[Vu{t') + F{t')]dt'. 



Using (4.1), (4.2) with s = and assumption (4.3) we get 



-1/2 



< 



\\f\\L2 + \\wl/^[Vu + F]\\L2L2 < 



If S is sufficiently small this implies both (4.4) and (4.5). 

To prove (4.7), (4.8) we proceed in a similar way using again (4.1) and (4.2): 

\\w-'/'Ku\\^.^. < ||A^/|U. + \\wl/'K{Vu + F)U.r.2. 

We shall need the following fairly elementary product estimate involving the angular 
derivative operator A^^ 

\\Ki9h)\\Liis2) < ||A^5||i2(s2)||A^/t||i2(s2) (4.9) 

which holds provided s > 1. This estimate can be proved e.g. by localizing the norm 
on the sphere via a finite partition of unity, and then applying in each coordinate 
patch a standard product estimate in the Sobolev space ^^^(R^), s > 1. 
Applying (4.9), and using assumption (4.6), we have 

\\w'J'K{Vu + F)U2L2 < S\\w-'/^Ku\\l2l2 + \\wl/^KFU.L2 

and the proof is concluded as above. □ 

We note the following consequence of (4.4): 

Corollary 4.3. Assume that the hermitian matrix V{x) satisfies, for 5 sufficiently 
small, C arbitrary and cr > 1 (with Wa-{r) = r(l + | logr|)°'j 

\vix)\ < -4n,, \^v{x)\ < ^ 



Wa{\x\)' ■ Wcr(\X\^ 

Then besides (4.4) we have the estimate for the derivatives of the flow 

ll-.-^/^Ve^*(^+^)/|L2,2 <||/|U.. 
// in addition we assume that, for some s > 1, 



\\KV{r-)\\ms2) 



< 



l|A:V\^(r.)||L^(s^) 



< 



c 



Wa\r) Wa{r) 

then we have the following estimate with angular regularity 
||u;;i/2vA^e^*(^+^)/|L2^2 <||A^/||^.. 



(4.10) 
(4.11) 
(4.12) 
(4.13) 
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Proof. Assume at first s = and let u — e'*^^+^''/. Each derivative uj = djU 
satisfies an equation like 

idtUj ~ Vuj + Vuj + VjU, Vj — djV, Uj{0, x) = fj = djf 

so we can represent it in the form 

Jo 

To the first term at the r.h.s. we can apply estimate (4.4) obtaining 

To handle the second term we use (4.5): 

< \\wI'^V,u\\l-2li < \\w„Vj\\l^\\w-'/^u\\i^2l2 

and again by (4.4) and by the assumption on WV we conclude the proof of (4.11). 

For the proof of (4.13) we apply to the equation for u the operator \D\ which 
commutes with V: 

idt{\D\u) = 'DilDlu) + \D\iVu) 
and we use estimates (4.7), (4.8), obtaining 

\\w-'/'\D\Ku\\^.^. < WKfWm + \U/'\D\K{Vu)U.^.. 

Now in the last term we commute \D\ with A and we notice that we can replace 
|D| by V obtaining an equivalent norm. This gives 

We can now apply the product estimate (4.9) and assumptions (4.12); proceeding 
as in the first part of the proof we finally obtain (4.13). □ 

By a similar perturbative argument, we obtain the endpoint Strichartz estimates 
for the Dirac equation with potential. Notice that in the version of this Theorem 
given in the Introduction (Theorem 1.2) we used an equivalent formulation in terms 
of the potential Vcr{x) = \x\2 \ log + (x)^'^"' . 

Theorem 4.4. Assume that the hermitian matrix V{x) satisfies, for 5 sufficiently 
small, C arbitrary and a > 1 (with Wa{r) = r{l + | logr|)'^j 

\n^)\ < , ,1/, 1,. ^ ivn^)i < , I,, . (4.14) 

{X)2+Wcr(\x\)^ (x)2+Wo-(|a;|)2 

Then the perturbed Dirac flow satisfies the endpoint Strichartz estimate 

If instead we make the following assumption (which implies (4.14) ): for some s > 1, 
\\KV{r-)U2^s.) < ^/ ||A^Vy(r.)||L2(s2) < ^ ^ ^^ , (4.16) 

then we have the endpoint estimate with angular regularity 

and the energy estimate with angular regularity 

IIA^e^^^^+^VlUrHi < WKIWh^ (4.18) 



'1/2 
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Proof. Consider first (4.15). Notice that V satisfies in particular the assumptions 
of Corollary 4.3. We can write 

with 

J = e'*^/, II = i[ e*(*-")^yMds. 

The term / is estimated directly using (2.36) with s = 0. On the other hand, 
applying (2.37) to the term // we get 

\\ii\\LiLrLi<m'^'-mvu)\\^.L.. 

Now we recall that the Riesz operators |Z3|~"'^V are bounded on weighted spaces 
with A2 weights, and {x)^ belongs to this class provided s < n/2 (see [21]). Thus 
we can continue the chain of inequalities as follows: 

= ||(x)^+|7^|-iV|i?|(F^.)|U.i. <||(x)^+V(F7.)|U.i. < A + B 

where 

A = \\{x)"^+{VV)u\\l2l2, B = \\{xY^+VVu\\l^l- 

Then we have 

A < \\{x)i+whv\\L^\\w^K\\L2L- < Wfh- 
by the assumptions on Vl^ and (4.4), while 

B < \\{x)^+wlv\\L^\\w^hu\\L2L2 < WfWm 

by (4.11). Summing up, we arrive at (4.15). 

The proof of (4.17) is similar. We estimate / using (2.3G). Applying (2.37) to 
the term // we get 

WKIIW^.^^^. < \\{xf2 + \D\A:,{Vu)h2L2. 

Then we commute \D\ with A^^, and we can replace the operator \D\ with V since 
the norm is equivalent; we arrive at 

\\KII\\lU-li < \\{x)i+Ki'^V)uU.L2 + ||(x)HA^y(Vu)|U.i.. 

Now we use the product estimate (4.9) and assumptions (4.16) to obtain 

< C\\w-'^^At,u\\L2L2 + S\\w-'/^KVuU2L2 

and recalling the smoothing estimates (4.7), (4.13) we conclude the proof of (4.17). 

It remains to prove (4.18). Consider first the free case V = 0. We have the 
conservation laws 

||e**^/|U^i2 = ll/lli., \\Ve''^f\\L^L2 = \\Vf\\L2 (4.20) 

which imply 

since ||2?/||i2 ~ Moreover, the operator A,^ introduced in (2.45) commutes 

with T), so that we have for all s > 

and switching back to the equivalent operator A^ as in the proof of Corollary 2.4 
we obtain 

\\Ke''^f\\L^m<\\Au.f\\m. (4.21) 
Consider now the case V ^ 0. We start from 

Uxyi-e^'-^fU^L^ < ll/IU. (4.22) 
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which is a consequence of (4.1) (we relaxed the weight). Taking the dual of (4.22) 
we get 



'Fit')dt' 



< 



{x)^+F\\l2l- 



L2 



which together with (4.20) gives 

e'^'-''^'^Fit')dt' 



< \\{x)^ + F\\l2l2. 

L°°L2 

Now a standard application of Christ-Kiselev' Lemma in the spirit of [13] (see 
also [4] for the case of Dirac equations) allows to replace the time integral with a 
truncated integral and we obtain 



s'(*-*')^F(t')dt' 



< \\(x)2+F\ 



L2L2 



L°-=L2 



(4.23) 



Recalling that the operator introduced in (2.45) commutes with 2?, and pro- 
ceeding as in the proof of Corollary 2.4 we obtain for all s > 







< 



-AtFl 



L2L2 



(4.24) 



and finally, applying \D\ which commutes both with V and Ai^, we have also 
ft 



Af 



s^(*-*')^F(t')d<' 



< \\(x)^+At\D\F\\r2 



L2L2 



(4.25) 



Now we use again the representation (4.19); by (4.21) and (4.25) we can write 

l|A:e»*(^+^)/||^^^. < WKfU^ + \\{x)^2+Kl\D\{Vu)U2^2 
and proceeding exactly as in the first part of the proof we arrive at (4.18). □ 

5. The nonlinear Dirac equation 

Theorem 4.4 contains all the necessary tools to prove global well posedness for 
the cubic nonlinear Dirac equation 



iut = Vu + Vu + Pi{u, u), m(0, x) = f{x). 
Our result is the following: 



(5.1) 



Theorem 5.1. Consider the perturbed Dirac system (5.1), where the 4x4 matrix 
valued potential V = V{\x\) is hermitian and satisfies assumptions (4.1G). Let 
P^{u^u) he a C^-valued homogeneous cubic polynomial. Then for any s > 1 there 
exists eg such that for all initial data satisfying 

WKfWm < eo (5.2) 
the Cauchy problem (5.1) admits a unique global solution u G CH^ Ci L^L°° with 

Proof. The proof is based on a fixed point argument in the space X defined by the 
norm 

\\u\\x ||A>||i2ic.i2 + ||A>|U=.^i. (5.3) 
Notice that in Theorem 4.4 we proved the estimate 

||e'*(^+^)/l|x<||A:,|Ui. (5.4) 
Define u = $(w) for v £ X as the solution of the linear problem 

iut = Vu + Vu + P{v,v), u{0,x)^f{x) (5.5) 
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and represent u as 

u = ^ e"(^+^)/ + i f e^(*-*')(^+^)p(«(t'),^)dt'. 

We recall now the product estimate 

\\K{9h)\\Li{s-) < |lA^g|li2(s2)|lA^/i|U2(s2) 
(see (4.9)). Then we have, by (5.4) 

\W\\x < WKfWm + r \\e'^'-'">'^P{v{t'),W))\\xdt' 



< WKfWm + I \\KP{v{t'),v{m\H^dt' EE WKfWm + \\KP{v,v)\\Lm^ 
By (4.9) we have 



whence 
and 



\\K{v')\\Lli^^) < l|A>||i2(s2) 



\\K{v')\\li<\\Kv\\li\\Kv\\Ul^^ 



\\K{v')hlLl < \\Kv\\l-li\\Kv\\%l^li < \Mx- (5.6) 
In a similar way, 

\\K^iv')\\Li{s-) < l|Af,Vi;|U2(s2)||A>||2,(g,) 

so that 

\\Kyiv')hi<\\K'^vh.\\Kv\\l^^. 

and 

\\K'^i^')\\LlLl < l|A^Vz;|U^^i2||A>|li2^^i2 < Ml^. (5.7) 
In conclusion, (5.6) and (5.7) imply 

\\KP{v,v)u.H^<\\vrx 

and the estimate for u = $(w) is 

||«||x^||<i>(t^)||^<||A^/||Hi + 111^11^. 
An analogous computation gives the estimate 

||$(w) - ^w)\\x < \\v - w\\x ■ {\\v\\x + \\w\\x f 
and an application of the contraction mapping theorem concludes the proof. □ 
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